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Numerical simulation of Fresnel diffraction with fast Fourier transform (FFT) is 
widely used in optics, especially computer holography. Fresnel diffraction with FFT 
cannot set different sampling rates between source and destination planes, while 



C/3 . 

Q , shifted-Fresnel diffraction can set different rates. However, an aliasing error may 

f/2 ' be incurred in shifted-Fresnel diffraction in a short propagation distance, and the 

^~>' aliasing conditions have not been investigated. In this paper, we investigate the 

aliasing conditions of shifted-Fresnel diffraction and improve its properties based on 

the conditions. 
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1. Introduction 

Numerical implementation of Fresnel diffraction [I] with fast Fourier transform (FFT) 
is widely used in optics, especially computer-generated holograms (CGH) and digital 
holography [21 [3]. CGH is a technique for generating a hologram pattern on a 
computer by simulating diffracted light from objects, and its applications include three- 
dimensional display and the design of optical elements. Digital holography reconstructs 
the light of an object including the amplitude and phase on a computer from a hologram 
captured by imaging devices such CCD cameras. 

In these fields, Fresnel diffraction and the angular spectrum method are used as 
key techniques. These diffraction calculation can be accelerated by FFT-implemented 
convolution; however it imposes the limitation that the sampling rate on a source plane 
is the same as that on a destination plane due to the property of FFT. If we can change 
the sampling rate, we will be able to extend the applications of diffraction calculations. 

In the angular spectrum method, some studies address different sampling rates 
on source and destination planes. Reference [1] proposed the scaled angular spectrum 
method using scaled-FFT, which is one version of FFT introducing a scaling parameter. 
This angular spectrum method requires four FFTs. Another one was proposed by us 
[5]. Our scaled angular spectrum method was implemented by non-uniform FFT. This 
angular spectrum method requires two FFTs and one interpolation. 

Fresnel diffraction to address the scale operation has already been realized by 
some studies. Double-step Fresnel diffraction executes the scale operation by double 
Fresnel diffraction calculations [6]. Shifted- Fresnel diffraction addresses the shift and 
scale operations by scaled-FFT [7]. Reference [8] proposed the same method as the 
shifted-Fresnel diffraction by Bluestein transform. Scaled-FFT and Bluestein transform 
are essentially the same. 

In holography, the scale operation of shifted-Fresnel diffraction [7j is used in several 
applications: for example, wavefront-recording method [9l [TOl [11] for accelerating CGH 
generation, digital holographic microscopy observable at any magnification [12], and 
zoomable holographic projection without a zoom lens [13]. Shifted-Fresnel diffraction 
yields new applications; however, the calculations may incur a serious aliasing error in 
a short propagation distance and the aliasing error has not been investigated. 

In this paper, we investigate the aliasing conditions of shifted-Fresnel diffraction 
and improve the problems incurred based on the conditions. The improved diffraction 
is referred as to aliasing-reduced shifted and scaled (ARSS)-Fresnel diffraction. In 
Section 2, we explain ARSS-Fresnel diffraction. In Section 3, we show a comparison 
of diffracted results between ARSS-Fresnel diffraction and shifted-Fresnel diffraction. 
Section 4 concludes this work. 
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2. Aliasing-reduced Fresnel diffraction with the scale and shift operations 

Fresnel diffraction is expressed by, 

U2{X2) = 7T / ui{xi)exp{ — {x2 -xi) )dxi (1) 

where k is the wave number, A is the wavelength of light, ui{xi) is the source plane, 
U2{x2) is the destination plane, and z is the propagation distance between the source 
and destination planes. 

We derive ARSS-Fresnel diffraction by the following relation: 

(x2 — sxi + o)^ = s{x2 — Xi)^ + (s^ — s)x\ + (1 — s)x\ + 2oX2 — 2soxi + 0^(2) 

where s is the scaling parameter, xi = prrii and X2 = pm2. mi and 7712 are integers with 
ranges of mi,m2 G [—N/2, N/2 — 1] where N is the number of sampling points, p is the 
sampling rate on the destination plane and is offset from the origin. The sampling 
rate on the source plane is determined by sp. 

We can obtain the following equation by substituting the relation into Eq.([l]): 

ui{xi) exp(— ((s^ - s)xl - 2soxi)) exp( )dxi 

AZ \Z 

(3) 

We can obtain ARSS-Fresnel diffraction by using the convolution theorem and 
introducing the band-limiting function to reduce aliasing: 



U2{X2) = C^T" 
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where Rect(-) is the rectangular function, x^ = pm^ ( nih G [—N/2, N/2 — 1]), and 
exp(i0„), exp(z0/i) and Cz are defined by, 

{s^ - s)xl - 2soxi . . 

exp(z0„) = exp(z7r ) (5) 

2 

exp{i(j)h) = exp{i7r—^) (6) 

AZ 

_ exp(z0,) _ exp{ikz + j^((l - s)x| + 2ox2 + o^)) 

^2 — ^T — ^T [') 
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Band-limitation function will be acceptable other than the rectangular function, for 
example, the circular function, window functions (e.g. Hanning and Hamming window 
functions) and so forth. It is straightforward way to extend to two-dimensional ARSS- 
Fresnel diffraction because Fresnel diffraction allows separation of variables. 

Figure [1] shows the light intensity distributions calculated by shifted-Fresnel 
diffraction [3 [8]. The calculation parameters are A^ = 1, 024, A = 633 nm and 2; = 0.1 m. 
The sampling rates on the source and destination planes indicate p\ and p2^ respectively. 
The scaling parameter s(= P1/P2) of Figs[I] (a), (b) and (c) are s = l{pi and p2 are 6 fi 
m), s = 6/4 and s = 6/10, respectively. 

As shown in the figure, shifted-Fresnel diffraction can magnify diffracted results by 
the scaling parameter; however, no countermeasure for aliasing has been devised. For 
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(a) 



(b) 



(c) 



Figure 1. Light intensity distributions calculated by shifted-Fresnel diffraction [71 |H]- 

(a) s = 1 (b) s = 6/4 (c) s = 6/10. 



instance, when changing z from 0.1 ni to 0.05 m and 0.03 m, the results are shown in 
Figl21 The upper images show the real parts and the bottom images show the light 
intensity of shifted-Fresnel diffraction. In the case of z = 0.1m, the results have good 
quality without aliasing. However, in the cases oi z = 0.05m and z = 0.03m, the results 
are extremely contaminated by aliasing. 



z=0.1m 



Real part 



Intensity 




Figure 2. Real part and intensity of shifted-Fresnel diffraction at the propagation 
distance of 0.1 m, 0.05 m and 0.03 m. 



2.1. Aliasing conditions 

Although ARSS-Fresnel diffraction is derived in a different way to shifted-Fresnel 
diffraction, ARSS-Fresnel diffraction of Eq.(|l]) is essentially the same as shifted- 
Fresnel diffraction, except for the introduction of the band-limitation function. In this 
subsection, we clarify the aliasing conditions of ARSS-Fresnel diffraction. Equation (jlj) 
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involves three chirp functions, namely exp(20c)5 exp(i0„) and exp(20/i). The cause of 
aliasing is the chirp functions. Figure [3] shows the real part of the chirp functions, under 
the calculation conditions of A^ = 1, 024, o = 0, A = 633 nm, p = lOyum and s = 6/10. 
In the case of 2; = 0.1 m, the aliasing is not incurred in the chirp functions; however, in 
the cases of 2; = 0.03 m, exp{i(j)c) and exp{i(j)h) are contaminated by aliasing. 



z=0.1m 



z=0.03m 




Figure 3. Real part of the chirp functions at the propagation distance of z = 0.1 m 
and z = 0.03 m. 



Therefore, we need to band-limit these functions. The aliasing-free areas of the 
chirp functions can be calculated by [1], 
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From Eqs.fl8|)-f lT0|) . the aliasing-free areas in pixel units are as follows: 
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According to Eq. (fT3|) . the band-limited chirp function of exp{i(j)h) sX z = 0.03m is 
shown in Fig|31 
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Changing z between 0.01 m and 0.25 m, FigJH plots the ahasing-free areas in pixel 
units of the chirp functions under the calculation conditions of A^ = 1,024, o = 0, 
A = 633 nm, p = lO/xm and s = 6/10. In the figure, the green, red, and light blue 
plots indicate the aliasing-free areas of exp(i0/i), exp{i(f)u) and exp(20c), respectively. 
The blue line indicates the half area {N/2) of the source and destination planes. The 
aliasing is incurred when each plot falls below the blue line. 

In FigJl](a), the effective distance without aliasing is over about 0.07 m because the 
chirp functions exp(i0c) and exp(i0u) do not incur aliasing over about 0.07 m, and we 
need to band-limit exp(i0/i) in about 300 pixels {x^ax ~ SOOp in Eq.(j4])). 

In FigHl^b), the chirp functions exp{i(j)u) and exp{i(f)h) for p = 20/ini and s = 6/20 
do not incur aliasing over about 0.18 m, while exp(z0c) still incurs aliasing at about 0.18 
m. If we band-limit exp(20c) in the aliasing- free area of about 200 pixels, no aliasing 
is incurred; however, the complex amphtude resulting in Eq.(j4l) outside the aliasing- 
free area was eliminated because zero area of exp(i0c) outside the aliasing-free area is 
multiplied with the result of inverse FFT of Eq. (jlj) . 

Therefore, for complex amplitude, we need to satisfy the following relation to avoid 
aliasing: 

|mi|,|m2| > iV/2 (14) 

For light intensity, we need to satisfy the following relation to avoid aliasing: 

|mi|>A^/2 (15) 

because the exp(z7r0c) can be ignored in the light intensity. The aliasing condition of 
the light intensity 1^2(2^2)! is mitigated compared to that of complex amphtude. 
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Figure 4. Aliasing-free area in pixel units of the chirp functions under the calculation 
condition of A = 633 nm. (a) s = 6/10 (p ^ XQ^iva) (b) s = 6/20 (p = 20^ni). 



In FigJSt^a), the effective distance without aliasing in the scaling parameter of 
s = 6/4 {p = 4/xm) is over about 0.02 m because the chirp functions exp(z0c) and 
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exp(i0tj) do not incur aliasing over about 0.02 m, and we need to band-limit exp(i0/i) in 
the aliasing- free area of about 250 pixels {xmax ~ 250p in Eq.(jl])). In Figj5]J^b), the chirp 
functions exp(27r0c) and exp(z7r0ft) in the scaling parameter of s = 6/2 {p = 2/im) do not 
incur aliasing over 0.02 m, while exp(i0u) incurs aliasing at the distance. If we band- 
limit exp{i(j)u) in the aliasing- free area of about 250 pixels, the aliasing is not incurred; 
however, the complex amplitude and intensity of Eq.(|l]) was eliminated outside the 250 
pixels. 
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Figure 5. Aliasing- free area in pixel units of the chirp functions under the calculation 
condition of A = 633 nni. (a) s ~ 6/4 {p ~ 4/im) (b) s = 6/2 {p = 2fiYn). 



3. Results 



Figure |6] shows the intensity distributions calculated by ARSS-Fresnel diffraction and 
shifted- Fresnel diffraction [71 E], under the calculation conditions of A^ = 1, 024, A = 633 
nm, p = 10/im, and s = 6/10. We change the propagation distance from z = 0.1 m to 
0.07 m and 0.04 m. As we can see, the results of ARSS-Fresnel diffraction do not incur 
aliasing, while those of shifted-Fresnel diffraction incur aliasing at z = 0.07 m and 0.04 
m. 

Figure [7] shows the intensity distributions calculated by ARSS-Fresnel diffraction 
and shifted-Fresnel diffraction, under the calculation conditions of A = 633 nm, p = 4/.im, 
and s = 6/4. As we can see, the results of ARSS-Fresnel diffraction do not incur aliasing, 
while those of shifted-Fresnel diffraction incur aliasing at 0.04 m. From Figsj6] and EJ 
ARSS-Fresnel diffraction can be applied to a wider propagation distance than shifted- 
Fresnel diffraction. 
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Figure 6. Intensity distributions calculated by ARSS-Fresnel diffraction and shiftcd- 
Frcsncl diffraction. The scaling parameter is s = 6/10. 
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Figure 7. Intensity distributions calculated by ARSS-Fresnel diffraction and shifted- 
Fresnel diffraction. The scaling parameter is s = 6/4. 



4. Conclusion 



We clarified the aliasing conditions of shifted-Fresnel diffraction and improved the 
diffraction, which was named ARSS-Fresnel diffraction. ARSS-Fresnel diffraction is 
useful for CGH calculation and digital holography because of the scale and shift property. 
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